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We discuss the effect of certain types of static disorder, like that induced by curvature or topo- 
logical defects, on the quantum correction to the conductivity in graphene. We find that when the 
intervalley scattering time is long or comparable to T0, these defects can induce an effective time 
reversal symmetry breaking of the hamiltonian associated to each one of the two valleys in graphene. 
The phenomenon suppresses the magnitude of the quantum correction to the conductivity and may 
result in the complete absence of a low field magnetoresistance, as recently found experimentally. 
Our work shows that a quantitative description of weak localization in graphene must include the 
analysis of new regimes, not present in conventional two dimensional electron gases. 
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Introduction. Graphene provides a two-dimensional 
electron system that is distinctly different from the two- 
dimensional electron gases (2DEGs) hosted in common 
semiconducting heterostrutures. The low energy elec- 
tronic states of graphene can be described by two sets of 
two-dimensional spinors associated to two independent 
points at the corner of the Brillouin zone (K and K'). 
In the absence of short range potentials the states asso- 
ciated to the different points are not coupled and, in the 
long wavelength limit, they satisfy the two dimensional 
Dirac equation. Theoretical calculations based on these 
assumptions [lj reproduce the unconventional quantiza- 
tion of the quantum Hall conductance observed experi- 
mentally (see alsoQ). It is expected that, beside 
the quantum Hall effect, other phenomena should mani- 
fest unusual characteristics as compared to more conven- 
tional 2DEGs. 

It appears that one of these phenomena is the weak- 
localization correction to the conductivity 0, 0, 0, @- I n 
the presence of time-reversal symmetry, the suppression 
of weak-localization at low magnetic fields produces a 
negative magnetoresistance, ubiquitous in metallic con- 
ductors. The magnitude of the weak-localization cor- 
rection is of the order of e 2 /h and it is determined by 
two characteristic time scales: the phase coherence time 
T0 and the elastic scattering time r. In the graphene 
samples that exhibit very clear QHE, however, no low- 
field magnetoresistance is observed. In these samples, 
nevertheless, the observation of high-order quantum Hall 
plateau indicates that phase coherent propagation of elec- 
trons occurs at least on a distance of several hundreds 
nanometers, which corresponds to the estimated elastic 
mean-free path, so that a magnetoresistance originating 
from the suppression of weak-localization is expected. 

It has been realized by Suzuura and Ando 9] that the 
quantum correction to the conductivity in graphene can 
differ from what is observed in conventional 2DEGs (see 



also|lfj|). This is due to the pseudo-spin associated to 
the solutions of the Dirac equation, that in conjunction 
with the nature of elastic scattering in graphene may 
change the sign of the localization correction, and turn 
weak-localization into weak-antilocalization (even in the 
absence of spin-orbit interaction) . The phenomenon cru- 
cially depends on the inter- valley scattering time t; v , i.e. 
the characteristic times that it takes for a charge car- 
rier to be scattered from one to the other iv-point in 
graphene. Specifically, if t; v ^> weak antilocaliza- 
tion is expected, resulting in a positive magnetoresistance 
at low field; if t 1v -C t^, conventional weak localization 
should be observed. Even though these conclusions are 
at odds with experimental findings, they are important 
as they clearly illustrate the relevance of inter- valley scat- 
tering. 

Here, we show that not only weak localization and 
antilocalization can appear in graphene depending on 
the nature of elastic scattering, but also that the quan- 
tum correction to the conductivity can be entirely sup- 
pressed due to time-independent potential slowly varying 
in space. These static potentials can result in the effective 
breaking of time reversal symmetry of electronic states 
around each if -point when rj„ > t^. They originate from 
defects that are realistically present in graphene, such as 
long-range distortions induced by topological lattice de- 
fects (disclinations and dislocations), non-planarity of the 
graphene layers, and slowly varying random electrostatic 
potentials that break the symmetry between the two sub- 
lattices of graphene. We estimate the effect of each one 
of these defects in terms of a characteristic time, which 
acts as a cut off for the time-reversed trajectories of elec- 
trons responsible for weak- localization phenomena. If the 
characteristic time associated to one of these mechanisms 
is much shorter than r^, weak- localization is largely sup- 
pressed. This explains the absence of weak-localization 
in the first magneto-transport experiments in graphene. 
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We emphasize that the potentials that we consider, be- 
ing static, do not actually break time reversal symmetry 
in graphcne. However, in the presence of these potentials, 
time reversal symmetry connects electronic states asso- 
ciated to the two different if -points. If the electron dy- 
namics is such that electrons cannot be transferred from 
one to the other if -point within their phase coherence 
time (i.e., if r Jt , > r^), the contribution to interference 
due to time-reversed trajectories (responsible for weak- 
localization) is suppressed pi). 

The model. We use the continuum approximation 
to the band structure near the K and K' points 
of the Brillouin Zone of graphene. The wavefunc- 
tions can be expressed as a two component spinor, 
[^A,K,t(r), ^b,k,-[(^)], where A and B stand for the two 
sublattices of the honeycomb structure. Other spinors 
can be defined for the K' point and the down spin ori- 
entation. The curvature of the graphene sheet can be 
included by generalizing the ordinary derivative operator 
to the covariant derivative [lilll^j l , fo llowing the standard 
procedure used for spinor fields [1J|. These effects are 
described below by the spin connection operator S C urv 
which breaks the effective time reversal symmetry around 
each K point. Similarly, certain topological lattice de- 



fects imply the existence of rings with an odd number of 
Carbon atoms. A dislocation has a pentagon at its core, 
and a dislocation requires a pentagon-heptagon pair. An 
odd numbered ring of Carbon atoms implies that the 
two sublattices are interchanged when traversing a path 
which encloses it. This needs to be taken into account in 
the continuum description, irrespective of the smoothness 
of the induced distortion. It can be described by a non 
Abelian gauge potential Edef 0, E3 • It has been shown 
that its inclusion allows for an accurate description of the 
electronic spectrum in curved graphene sheets, such as 
fullerenesflil Il3| . or graphitic cones0. Note that this 
gauge field interchanges the two sublattices and the K 
and K' points. This fact, however, does not prevent the 
description of the system in terms of two sets of spinorial 
wavefunctions which obbey equivalent Dirac equations. 
In fact, using the transformation: 

*AKk,( f ) = Vb( f ) + ! *Mi s ( f ) 
*BK<k s ( f ) = V'L( f )-^Mk s ( f ) C 1 ) 

the model can be reduced to two independent hamiltoni- 
ans: 



A(f) v F (id x + d y ) + A x + iA 

v F (id x - d y ) + k x - iky + S curv + S def -A(r) 



^ = / A(f) v F (id x -dy)+A x -iA y + 'E CUIV -'E de { \ 

K v F (id x +dy)+k x + iky + 'E cul . v -'E dc[ -A(f) / 



where A is the ordinary vector potential. A(f) = 
Va (?) — Vq (?) measures the difference of the electrostatic 
potential in the two sublattices and originates from the 
environment around the graphene layer (e.g., charges lo- 
cated at random position in the substrate supporting 
graphene). 

A description of quantum interference in terms states 
associated to two decoupled K-points separately is cor- 
rect if t; v 3> t^. In this regime, the defects described 
by three terms Edcf , S curv and A(f ) break the effective 
time reversal symmetry around the K and K' points, 
which can be written as the complex conjugation oper- 
ator times cry^l. We will treat separately the effect of 
the different defects by analyzing the dynamical phase 
that quasi-classical trajectories in real space acquire in 
their presence. We will confine ourselves to the qualita- 
tive analysis of this problem, but we note that a formal 
description of weak -(anti)localization in terms of quasi- 
classical trajectories |g, M, lal can be made equivalent to 
diagrammatic approaches [17] . 




FIG. 1: Geometry of the scattering process analyzed in the 
text. An impurity potential is placed at the center of a spher- 
ical bump in a flat sheet. 



Curvature effects. Scattering is modified in a curved 
surface. We analyze the geometry sketched in Fig.pQ. 
An incoming plane wave moves into a spherical cap 
parametrized by the angle 9. The wave is scattered by 
a defect within the bump. We place the defect at the 
center of the region, in order to simplify the calculation. 
In flat space, the incoming and outgoing waves can be 
written as: 

* k - ( X ) ^ s ( X ) ^ 

where <f>^ determines the direction of k. The scattering 
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and Sdef in eq.© becomes: 




FIG. 2: Wavefunctions within the spherical cap shown in 
Fig.pP for 1 — 20. These two solutions, in flat space, are re- 
lated by the effective time reversal symmetry defined in each 
valley. 



amplitude due to a local potential, V = VS(r) is A = 

<k|i>|k')ocVqi + cos(<^-<^,)]- 

We calculate the scattering amplitude in the geom- 
etry shown in Fig.pQ by matching the plane waves in 
eq.@ to solutions inside the spherical cap. The wave- 
functions can be computed analytically for energies e; = 
w/Ry/l(l + 11 — 2|l2|. The quasiclassical limit corre- 
sponds to 1 3> 1. The main difference between the so- 
lutions inside the cap and those in a plane, eq.®, is that 
the values of the modulii of the two components of the 
Dirac spinors are not equal. Because of the radial sym- 
metry of the problem, we only need the solutions inside 
the cap with m — 0. The relevant solutions, for a given 
valley, can be written as: 



(4) 



Matching these solutions to those in eq.(@ outside the 
cap, the scattering amplitude can be written as: 

Acxy[|^ 1 4 (0 o ) 2 +^|(0 o ) 2 + 2^(0 o )V|(0o)cos(0 £ -^,)] 

(5) 

where 8 is the angle which defines the shape of the cap. 
The deviation from scattering in a plane is determined by 
the difference ip\(9o) — ^b^o), which is maximum when 
^i(0o), ip%( S o) « 0, as shown in Fig.0. 

The first zeros of the functions ip\(0o),tp%(0o) ~ lie 
at 6q ~ l - ■ Hence, we expect deviations from antilocal- 
ization for Fermi energies such that 9q ~ epR/vp, that is, 
kpR ~ 1. In a sample with height fluctuations of order 
h oner a length I, we have R^ 1 ~ h 2 /l, so that antilocal- 
ization effects will be suppressed for (kph?)/l ~ 1. 

Gauge field induced by lattice defects. As mentioned 
above, a local rotation of the axes leads to the existence 
of a gauge field, £def- In terms of the lattice strain, 
u x (r), u y (r), this rotation is given bv[is|: 



9(r) 



s dcf = -ve(r) 



(7) 



We consider first the effect of gauge field induced by a 
finite distribution of dislocations, with random Burgess 
vectors hi, and average distance d, on a quasiparticle 
moving around a closed loop of length I. 

At sufficiently long distances from an edge dislocation 
with Burgess vector b, eq.|@ gives 0(r) oc br/|r| 2 . The 
local rotations average to zero when the loop encloses 
completely the dislocations. A finite contribution is ob- 
tained from the dislocations whose core is near the tra- 
jectory of the particle. The quasi-classical width of the 
trajectory of the particle is of order k^ 1 , where fcp is 
the Fermi wavevector. Hence, the number of dislocations 
which contribute to change the phase around the loop is 
of order (k^ 1 d)/l 2 1 where I is the length of the loop and 
d is the average distance between dislocations. 

To estimate the effect of the rotation on the phase 
of the trajectories traversing the loop, we assume that 
the Burgess vector of the dislocations is distributed ran- 
domly. Using the central limit theorem, we find that 
the rotation associated to the loop becomes non neg- 
ligible when [(fcp 1 d)/l 2 ] 2 3> 1. This rotation leads to 
a phase of order it, which dephases randomly trajecto- 
ries traversed clock- and anticlockwise, which suppresses 
the quantum interference correction to the conductivity. 
From the above inequality, the time scale at which these 
effects are relevant is: 



Vp 



gauge 



k^d 



(8) 



Potential gradients. The potential A(f) in eq.Q rep- 
resents the difference in potentials at the two sites of 
the unit cell. Physically, this potential originates from 
char ges located in the substrate supporting the graphene 
laver|l9j. An asymmetry between the two sublattices 
arises from slowly varying potentials, that can be writ- 
ten as the sum of a smooth term which is the same within 
each unit cell, and a small contribution which breaks the 
equivalence of the two sites in the unit cell. The second 
part can be written as A(r) « cVV(r), We assume this 
potential to be slowly varying, so that: 



V(q) 



V 




|q| < q c 
|q| > q c 



(9) 



where q c is a momentum cutoff such that q c <c a -1 , where 
a is the lattice constant. We are interested in momentum 
transfers |q| ~ kp, where kp <C a -1 for typical dopings. 
Eq.ljnj implies that: 



(6) 



A(q) 



qcVb 




|q| < 9c 
|q| > q c 



(10) 
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The elastic scattering time can be obtained from Vdisorder 
using Fermi's golden rule: 



V 2 v(e F ) 



clasli< 



(ii) 



where f(ep) = |eF|/(7r« F ) is the density of states at the 
Fermi energy. 

The potential A(f) breaks the effective time reversal 
symmetry around each K point. We can estimate the 
inverse time at which this operator induces a significant 
change in the wavefuction of the spinor applying again 
Fermi's golden rule to the symmetry breaking potential: 



'grad 



ftp 



clastic nc l oc 



(12) 



where n c ioc ~ (kpa) 2 is the number of carriers per unit 
cell (note that the effect can be larger close to the sample 
edges where the disorder is also larger) . In typical exper- 
iments, n c \ ~ 10~ 4 — 10~ 2 . Thus, there exists a range 
of energies or temperatures, r~* d < T < tJ 1 , t^ 1 < T, 
where the effective time reversal symmetry around each 
K is not broken by A(r). 

Conclusions. Our analysis illustrates that the behav- 
ior of the quantum correction to the conductivity in 
graphene is much richer than what was anticipated in 
the work of Suzuura and Ando || . We conclude that the 
inter-valley scattering time Ti V and the phase coherence 
time Trf, are not the only important time scales. An ad- 
ditional time scale describing the effective time reversal 
symmetry breaking is present, which can cause a com- 
plete suppression of weak (anti)localization. This time 
scale depends on specific defects present in the graphene 
samples, which leads to the prediction that large differ- 
ences in the quantum correction to the conductivity mea- 
sured on different samples should be expected. This is 
striking, since in all metallic conductors studied in the 
past weak (antilocalization manifests as a robust and very 
reproducible phenomenon. 

Furthermore, the physical understanding provided by 
our analysis, which does not rely on any detailed assump- 
tion, enables us to draw additional conclusions. For in- 
stance, we expect that it will be easier to observe weak- 
localization in narrow graphene samples, since scattering 
at the edges couple states at the two different K-points. 
We also expect that an effective time reversal symmetry 
breaking similar to the one discussed here should occur 
more in general, in systems with two or more degenerate 
valleys with topological defects similar to those consid- 
ered here. This is the case, for instance, for a graphene 
bilaver[20l Elj ]. Note, however, that contrary to indi- 
vidual graphene layers, in graphene bilayers the relative 
phases of the two components in the spinor of the momen- 
tum eigenfunctions are twice those in graphene. There- 
fore, in sufficiently clean samples (where defects are not 
sufficient to induce an effective time reversal symmetry 



breaking) the usual negative magnetoresistance is ex- 
pected, and there should be no weak antilocalization a 
la Suzuura- Ando. 
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Note added. Since this work was posted, three related 
works have appearedjU |24| dealing with the same 
topic. Ref.|22j presents experimental results consistent 
with our work, as well as a theoretical explanation along 
similar lines. The analysis of weak localization effects 
inj^ is also consistent with our results. An alternative 
explanation is discused in|23l| . 
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